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We present wide-field, high-resolution imaging observations in ^^CO 3^2 and H2 1- 
S(l) towards a ~1 square degree region of NGC 2264. We identify 46 H2 emission objects, 
of which 35 are new discoveries. We characterize several cores as protostellar, reducing the 
previously observed ratio of prestellar/protostellar cores in the NGC 2264 clusters. The length 
of H2 jets increases the previously reported spatial extent of the clusters. In each cluster, 
<0.5% of cloud material has been perturbed by outflow activity. A principal component anal- 
ysis of the ^^CO data suggests that turbulence is driven on scales >2.6 pc, which is larger 
than the extent of the outflows. We obtain an exponent a=0J4 for the size-linewidth rela- 
tion, possibly due to the high surface density of NGC 2264. In this very active, mixed-mass 
star forming region, our observations suggest that protostellar outflow activity is not injecting 
energy and momentum on a large enough scale to be the dominant source of turbulence. 

Key words: ISM: jets and outflows - infrared: ISM - submillimetre - stars: formation 



1 INTRODUCTION 

Studies of young clusters are essential for understanding star for- 
mation, since most stars, and especially massive stars, are known 
to form in clusters (Zinnecker & Yorke 2007). NGC 2264 is an at- 
tractive target for studying this mode of star formation, since it is 
nearby, has relatively low foreground extinction and there is a large 
population of pre-main sequence stars, whose age spread is evi- 
dence of sequential star formation ( [Adams, Strom & Strom|1983]). 
The region contains a wealth of submm cores ([Ward-Thompson et 



aTpOOO^, pre-main sequence stars and protostellar sources ( [Young 



et al.[[2006{[Teixeira et al.[[2006|), many of which have formed in 



clusters ( [Wolf-Chase et al.|2003| [Furesz et al.'2006'; l Young et al!] 
|2006}. At 760-900 pc distance (see [Baxter et al. 2009} and refer- 
ences therein for discussion), the region contains over 300 Class 
II objects, and 30 Class I objects ( [Dahm & Simon|2005|). S everal 
young Class obje cts have been identified (iTeixeira et aL[[2006[ 
[Furesz et all2006[ [Peretto, Andre & Belloche"2006'). Kinematic 
motions include large-scale infall motions ( Peretto et al. 2006, 
[Williams & Garland|2002[), and large, fast ou tflows, including the 
well-known NGC 2264 G ^Lada & Fich|T996l[Hedden et al.[2006[ ). 
The wide field optical images of [Reipurth et al. ( [2004[ ) reveal a 
number of HH objects, several of which trace giant, parsec-scale 
outflows. Published infrared images of the NGC 2264 region cover 
only very small areas; maps of the NGC 2264 G molecular out- 
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flow (MHO 1358/135^ have been presented by [Davis & Eisloffel 
( [1995| , while more recently a number of H2 features (MHO 1349- 
1359) have been identified around NGC 2264 C by [Wang et"aLl 
( .2002J . 

Outflows have a significant impact on material within the par- 
ent molecular cloud, injecting energy and momentum at large dis- 
tances. The outflows extend typically over 0.1-1 parsec, and the 
combined sub-arcsecond resolution H2 images with deep spectral 
imaging in ^^CO, result in a more complete census of star forma- 
tion activity in clustered regions. Outflow kinetic energies can be 
comparable to, or even larger than, turbulent and gravitational en- 
ergies of natal clouds, and could have a significant disruptive effect 
( [McKee & Ostriker|2007| ). In order to find evidence for, and quan- 
tify the disruptive effects of molecular outflows, it is crucial to map 
the entire cloud. 

Fig. [1] shows the Digitized Sky survey (DSS) optical image 
of NGC 2264, with the main star formation regions and nebulae 
labelled, obtained with Starlink Gaia using SkycarJA detailed re- 
view of this region has been published by'Dahm[( [2008[ ). We present 
wide field, high resolution mapping observations of NGC 2264. 
Through CO observations we can provide a complete census of 
outflow activity in the region, investigating any dynamic impact 



^ The MHO (Molecular Hydrogen emission-line Object) catalogue is 
maintained at http://www.jach.hawaii.edu/UKIRT/MHCat/ 
^ http ://archive.eso.org/cms/tools^documentation/skycat 
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Figure 1. DSS optical image of NGC2264 (the Christmas Tree Neb- 
ula) showing the massive O star SMon, the two star-forming clusters 
NGC2264 C (also known as Allen's Source, and IRSl) and NGC2264 D 
(also known as IRS2), and the two nebulae within the cloud, the Cone Neb- 
ula and the Fox Fur Nebula. 
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Figure 2. Integrated CO intensity map of NGC 2264, in the velocity range 
-20 to 30 kms"\ with contours at 15,30,60,90,120,150 Kkms"^ A dia- 
mond marks the position of SMon. Crosses mark dense cores identified 
through SCUBA observations ( |Di Francesco et al^|2008t . Labels mark 
the well-known protostar NGC 2264 G, and the two star-forming clusters 
NGC 2264 C and NGC 2264 D. 



the outflows are having on the parent molecular cloud through en- 
ergy injection. Using H2 narrow line imaging, we can detect the 
youngest flows, from emission arising in gas shocked by the im- 
pact of protostellar outflows, and measure the extent of these flows 
by tracing the detected emission knots and HH objects. 



2 DESCRIPTION OF OBSERVATIONS 
2.1 Spectral Line Observations 

The spectral-line observations were taken at the JCMT (James 
Clerk Maxwell Telescope) using HARP (Heterodyne Array Re- 
ceiver Programme, |Buckle et al.|2009] ). Observations of CO 3^2 
at 345.796 GHz were taken with a 0.488MHz (0.423 km s"^) spec- 
tral resolution and gridded with 6" spatial pixels. The beam size is 
14'' at 345 GHz. The map covers an area of 1 square degree, and 
observations totalling 12.9 hours were taken during August and Oc- 
tober 2007, in position-switch raster mode. The observations were 
taken in Band 3 weather, with zenith opacity values between 0.13 
and 0.19, and average system temperatures ~550 K. Averaged over 
emission free regions of the map, a Icr rms of 0.32 K per pixel in 
1.0 kms~^ spectral channels was obtained. 

All intensities reported here are in units of TJ^, the an- 
tenna temperature corrected for sky and telescopes losses ( |Kut-| 
|ner & Ulich,1981| ). The main beam brightness temperature, Tmb = 
T^/y/mb. with a main beam efficiency, z/mb = 0.66. Frequent point- 
ing and focus observations were carried out, and calibration obser- 
vations using N2071IR. The scatter in these observations suggest a 
calibration accuracy of ~20%. 

The data were reduced and analyzed using the Starlink project 
software, in particular SMURF ( |Jenness et aL|2008| ) and KAPPA 
HCurrie et al.|2008| routines. 



2.2 H2 Imaging Observations 

The infrared imaging observations were taken at UKIRT (the 
United Kingdom Infrared Telescope) using the near-IR wide-field 
camera WFCAM ( Casali et al. 2007 ). The observations cover an 
area on sky of 0.75 sq. degrees (a WFCAM 'tile'), with a pixel 
size of 0.2''. Images were obtained through broad-band K (A = 
2.20 |j,m, 6 A = 0.34 \im) and narrow-band H2 (A = 2.121 \xm, 
6A = 0.021 \im) filters. Exposure times of 5 sec and 40 sec were 
used; thus, the total per-pixel integration times were 100 sec and 
800 sec in K and H2, respectively. 

Data were initially reduced by the Cambridge Astronomical 
Survey Unit (CASU), and distributed through the Wide Field As- 
tronomy Unit (WFAU) archive. Further data reduction and analysis 
was carried out using Starlink KAPPA routines ( ,Davis et al.,2009| l. 



2.3 Data Summary 

The two datasets presented here probe the most active regions of 
current star formation. CO 3^2 emission traces dense, warm gas 
typical in regions of star formation, with temperatures 10-50 K and 
gas densities of lO'^-lO^ cm~^. The higher critical density of the 
3^2 line over lower excitation lines makes this transition a use- 
ful tool for tracing dense, more collimated flows from the youngest 
sources. The rovibrational H2 1-0 S(l) line at 2.122 ^tm arises 
in regions with temperatures ~2000-4000 K, and densities 10^- 
10^ cm~^. Emission from this line indicates the presence of deeply 
embedded driving sources, and traces the youngest outflows. 

Fig. |2] shows the integrated intensity map of CO 3^2 emis- 
sion, in the velocity range -20 to 30 km s~^ (note all velocities are 
LSR unless otherwise stated). The brightest integrated emission 
(246 Kkm s"^) is towards the cluster NGC 2264 C, while the bright- 
est peak intensity (28 K) is towards the head of the Cone Nebula. 
Filaments can be seen extending south from NGC 2264 C, east and 
west from NGC 2264 D and north-west from the S Mon. S Mon is a 
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massive (17.8 M©) young star, of spectral type 07V ( [Herrero et al.| 
|1992| ), and the radiation pressure and wind from this star is hkely 
to be heating and dynamically affecting nearby regions, such as the 
protostellar outflow source NGC 2264 G ( Teixeira et al.|2008[|Ladal 
|& FichI 19961 ). The two outflow lobes of NGC 2264 G can be seen 
to the north-east of S Mon. West of S Mon, and extending south- 
wards, is a region of very weak emission, which we refer to as a 
cavity, surrounded by a ring or bubble of emission containing sev- 
eral SCUBA dust cores. 

Dust pillars and cometary clouds, harbouring YSOs, are 
prevalent in Spitzer imaging of massive star forming regions ( |Smith| 
|et a l. 2010), where the emission is dominated by strong PAH fea- 
tures. The Cone Nebula resembles one of these cometary regions in 
CO emission. The structure of this emission suggests that molecu- 
lar gas is being ablated, possibly by the action of the NGC 2264 C 
cluster, exposing star forming cores. Red- and blue-shifted ^^CO 
emission and H2 emission at the head of the Cone indicates on- 
going star formation activity. 

Fig. [3] shows a colour composite image of the H2 1-0 S(l) 
and K-band images of NGC 226^ Ellipses mark the Molecular 
Hydrogen emission-line Objects (MHOs) we have detected. Inset 
images show zoomed versions of these data towards two regions, 
the cluster NGC 2264 D and an isolated outflow to the south-east. 
H2 line emission appears as turquoise in these images, which are 
shown in detail in Appendix Figs. [TT] to [14] We detect 46 MHOs, 
of which 35 are new discoveries, and described in detail in the 
Appendix. We detect regions of multiple jets and bow shocks sur- 



rounding NGC 2264 C (extending further than seen by |Wang et 
200?), NGC 2264 D and NGC 2264 G (as seen by [Reipurth et 



2004). There are additional regions of jet activity to the east of 
NGC 2264 D, and in the south-east corner of the map. Towards the 
Fox Fur Nebula and the Cone Nebula, we see much larger- scale 
arcs and filaments in H2 emission. 

Representative CO spectra, and the intensity- weighted veloc- 
ity map are shown in Fig.|4] highlighting the varying and complex 
velocity structure of this region. To the north, emission from the 
Fox Fur Nebula has narrow linewidths, and velocities >9 kms~^ 
The star forming clusters, and NGC 2264 G have extended line 
wings due to outflow activity. The clustered regions of star forma- 
tion are kinematically distinct, with different cloud velocities, as 
has been previously noted (Maur y et al.|2009[|Peretto et al.|2006 l. 
The Cone region, and NGC 2264 C, have velocities ~7 kms~ , 
while NGC 2264 D has velocities ~5 km s~^ Several regions show 
multiple velocity components, including the cavity south-east of 
S Mon, the Cone Nebula and NGC 2264 D. 

ICrutcher, Hartkopf & Gig"uerel ( [T978l ) and |Furesz et"aL| ( [20061 ) 
identify distinct regions across NGC 2264 which are broadly in 
agreement with the velocity gradient and extent of the red, green 
and blue velocity regions in Fig. |4] respectively. Sung, Stauffer & 
|Bessel| ( |2009| identifies and assigns ages to distinct star forming 
clusters using Spitzer data. The region surrounding S Mon, includ- 
ing NGC 2264 G is aged at 3.1 Myr. The Cone Nebula is at a 
younger age, due to the number of embedded sources. [Peretto, HelT] 
[nebelle & And re ( 2007 ) calculates very young ages for the clusters 
NGC 2264 C and NGC 2264 D of ~0. 1 Myr. We look at the regions 
of active star formation in more detail in the following section. 
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available from 



|http://www.jach.hawaii.edu/UKIRT/TAP/| 



3 ACTIVE STAR FORMING REGIONS 

The wide-field ^^CO 3 ^ 2 and H2 1-0 S(l) data can be used to 
provide a comprehensive quantitative analysis of the star formation 
activity across NGC 2264. Outflow and jet properties are expected 
to change as proto stars evolve from the early accretion phase. The 
youngest objects, which are still accreting most of their final mass, 
are expected to have more powerful outflows, as measured from 
the CO momentum flux or jet luminosity ([Caratti o Garatti et aL| 
[2006 , Arce & Sargent 2006)|Bontemps et al. 1996 ). In the following 
analysis, we use CO 3 ^ 2 to estimate the total outflowing mass, 
and the energy and momentum injected into the ambient cloud. We 
use H2 1-0 S(l) to estimate the total H2 luminosity in the jets, and 
place constraints on the impact of star formation activity on the 
ambient cloud. 

For the emission in the outflow wings, we assume optically 
thin conditions, and calculate the total mass and energetics of 
molecular gas using emission from the ^^CO 3^2 transition 
following Garden et al.| ( [r991) ; Buckle et al.| ( [20T0l ). [Maury et al.] 
(2009 ) examined the opacity in the line wings of ^^CO using CO 
isotopologues towards the NGC 2264 C cluster, and found the emis- 
sion to be optically thin outside of the central ~10 km s~^ of the line 
profile. For the outflow contribution at low velocities, the optically 
thin assumption may not be correct, and the masses are then lower 
limits. Since we make no correction for the inclination angle of the 
outflows, the values for momentum and kinetic energy are lower 
limits. We use an excitation temperature of Tex = 20 K in order to 
aid comparisons with previously published work (e.g. |Lada & Fich| 
[T996J . The data are thresholded at the 3cr level. 

Due to the NGC 2264 velocity field (shown in Fig. [4|, care 
needs to be taken in assigning velocity limits for the high veloc- 
ity emission. We have taken the values from Peretto et al.| ( |2006] ), 
who describe N2H+ observations, giving systemic velocities of 
7.5 km s"^ for NGC 2264 C and 5.5 km s"^ for NGC 2264 D. For the 
high velocity material, we use the velocity intervals +5.0 kms~^ to 
±20 km s~^ either side of the systemic velocity for each cluster. For 
other regions, the sources are sufficiently isolated for an estimate 
of the systemic velocity to be made from ^^CO emission in areas 
adjacent to the molecular outflows. For NGC 2264 G, we use a sys- 
temic velocity of 5.0 km s~^ , with ranges up to ±30 km s~^ from the 
systemic velocity. Emission from ambient material covers the red 
lobe at low velocities, and so we use a lower velocity of 2 kms~^ 
for the red lobe, and -1 km s~^ for the blue lobe. For the protostar 
in the far south, IRAS2, we adopt a systemic velocity of 7.8 km s~^ , 
and velocity intervals of 2 to 10 km s~^ for the red- shifted emission, 
and -6 to -1 km s~^ for the blue-shifted emission. 

We have extracted H2 line emission fluxes from the H2 plus 
continuum data after carrying out aperture photometry with back- 
ground subtraction to calculate the calibrated integrated line flux. 
The H2 plus continuum data is used rather than the continuum- 
subtracted data, since it does not add to the calibration uncertain- 
ties, and there are no deep negative components from an imperfect 
subtraction routine. Calibration was checked on stars in the field 
with known K-band magnitude, and is accurate at the 10% level. 

The quadrupole transitions of H2 (known as O, Q and S 
branches) are the main cooling mechanism for shock-excited jets 
detected in the near-IR. The rovibrational transitions arising in the 
near-IR, including the H2 1-0 S(l) transition, can be used to es- 
timate the total molecular hydrogen luminosity (L//2) of the flow. 
Younger flows are expected to have more powerful outflows, and 
higher values of Lh2 • 

The total H2 luminosity (LH2) can be calculated from the H2 1- 
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Figure 3. Colour composite image of the H2 + continuum (blue and green) and K-band (red) data sets. Narrow-band H2 line emission is seen in turquoise. 
The contour shows integrated CO emission at 15 Kkms~^. Ellipses mark the positions of MHOs we have detected, described further in the Appendix. Insets 
show zoomed images of two regions of newly detected MHOs, NGC 2264 D and a region in the south-east. 



S(l) luminosity (Lh^ 12). with Lh^ ~ IOXL//2 12 ( [Caratti o Garatti et| 
|al.|2006| ), for typical H2 jet temperatures between 1500 and 2500 K. 
If emission from the jet arises in regions at higher temperatures, this 
method will underestimate Lh2 by factors up to 2.5. Lh2 will also 
be underestimated by up to an order of magnitude in regions of 
high extinction (Av=25). |Caratti o Garatti et al.| ( |2006] ) have used 
near-IR spectral line observations to calculate the temperature and 
extinction towards one of the flows in our dataset, NGC 2264 G. 
Their values of T=2 100-2800 K, and Ay =3-8 mag suggest Lh2 = 
10 X Lh2 12 is a reliable estimate of L//^ for this region. 



Table [T] shows the mass, momentum and kinetic energies cal- 
culated for the red- and blue-shifted CO 3^2 emission, and the 
H2 luminosities, towards each region. Individual MHOs included 
for each region are described in the Appendix. For the whole cloud, 
L//2 = l-l L©, although this is a lower limit, since we have not taken 
into account extinction. The masses are lower limits, since they are 
dependent upon temperature and opacity. If the high velocity emis- 
sion is excited at temperatures >20 K, or the assumption of opti- 
cally thin emission is not valid, then the masses will be increased. 



Wide-field imaging of NGC 2264 5 






-10 


-5 





5 


10 


15 20 25 


L2 
8 








ivJ 


\ 


NGC2264 D 


4 






/ 




V 







r^^^^^^-^^ 




^s^ 






-^^-UKT^.^-,,..^...^ 





-10 


-5 





5 


10 


15 20 


25 


25 








jl 




Cone 


- 


/ 20 














- 


^ 15 








f \ 






- 


^10 










^ \ 




- 


- 5 










\ 




- 





—— ^xu-~ 


—, 


— r^ 
















-10 -5 5 10 15 20 25 
Velocity (kms"M 



41m 6h40m 

Right ascension (J2000) 



Figure 4. Spectra from individual pixels, from locations shown in the intensity-weighted velocity image, 2 TtVi/Av 2 Tt, where Tt, vt, Av/ are the intensity, 
velocity and velocity channel width respectively for each velocity channel, which is used to determine the velocity field. The intensity-weighted velocity data 
are thresholded to the >3cr detection limit. The colour scale indicates the velocity field in units of km s~^ . 



Additionally, it is not possible to entirely separate emission from 
outflowing gas from that of the ambient cloud. 

Towards all of the active star forming regions, we find more 
mass, momentum and energy in the red- shifted material than in the 
blue-shifted material. 'Maur y et"ar] ( |2009| ) has examined the indi- 
vidual flows in NGC 2264 C using CO and isotopic data, and also 
finds the red- shifted material to contain more mass and momen- 
tum. Asymmetry in molecular outflows is relatively common in low 
mass protostars, although no mechanism has so far been proposed 
for the brightening of the red- shifted lobe. 



3.1 The NGC 2264 protoclusters 

NGC 2264 D is a more massive cluster than NGC 2264 C, with 
more protostellar sources, more mass and energy, and covering a 
larger area. Within the uncertainties, both clusters have similar H2 
luminosities. Fig. p] shows the red- and blue-shifted CO emission 
contoured over the H2 image towards NGC 2264 C. The jets and 
outflows in the cluster can be seen extending past the Cone Nebula 
to the south. To the south-west of IRSl is a red-shifted flow that 
contains the most high velocity gas, with line wings detectable out 
to ~30 km s~^ . We have described in detail the individual H2 objects 
and CO flows in the Appendix. 

[Peretto et al] ( 2006 ) identify 4 of the 12 NGC 2264 C cores 
as prestellar. As described in the Appendix, all four of these cores 
are spatially located within collimated red- and blue- shifted CO 
emission, and compact knots of H2 emission. Additionally, 3 of the 
cores are coincident with Spitzer 24 \xn\ sources. These data suggest 
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^ Emission from MH01349 is excluded, since it can't be separated from 
the saturated continuum towards Allen's source. 

Table 1. Outflow parameters for the NGC 2264 high velocity regions. 



that these 4 cores are young Class O/I protostars. |Sung et al~|p009| l 
classifies 17 24 |xm sources as protostellar in NGC 2264 C. 

[Peretto et al.| ( |2006| ) estimate the total gas mass in 
NGC 2264 C to be 1650 Mo, compared to the 3 Mo we estimate 
to be entrained in the outflow. Although there are multiple ener- 
getic outflows in this cluster, only ~0.2% of the molecular gas has 
been entrained by outflows, or aflfected by star formation activity. 
This is a lower limit, since it is dependent upon the opacity. 
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Figure 5. Red and blue shifted CO emission contoured over the H2 continuum- subtracted greyscale image towards the NGC 2264 C cluster. Contours are at 
3.9, 6.5, 9.1, 14.3, 19.5, 27.3, 37.7, 48.1 Kkms~^ Red-shifted emission extends from 12.2 to 30.9 kms~\ and the blue-shifted from -17.8 to 2.1 kms~^ The 
saturated source seen in the H2 emission is Allen's source jAllen|1972t , IRSl. 



Fig. |6] shows red- and blue-shifted CO emission contoured 
over the H2 image towards NGC 2264 D. The cluster extends fur- 
ther westwards than previously observed, with the discovery of new 
jets and molecular outflows spatially adjacent to this cluster. The 
red-shifted emission contoured at the south-eastern edge of the map 
is from the NGC 2264 C cluster. The total gas mass is 1310 M© in 
NGC 2264 D ( |Peretto et aL|2 006), so the molecular outflows, en- 
training 5 Mq, constitute a negligible fraction (~0.4%) of the clus- 
ter mass. Although this is a lower limit, this result implies that the 
large number of outflows are not aff'ecting the bulk of molecular 
material. 

I Sung et al. | ( |2009| ) classifies 44 sources as protostellar in 
NGC 2264 D. There are many overlapping CO outflows and H2 jets, 
which we are unable to separate near the sources, making source 
identification difficult. [Forbrich et al.| |2010) identifies 3 embed- 
ded YSOs as candidate massive YSOs using PAH emission. One 
of these, SSB 11829 is co-incident with a ^^CO bipolar molecular 
outflow, and MHO 1385, a bright knot and extended west-facing 
bow shock, indicating the source is protostellar. A more detailed 
description of the individual H2 objects and CO flows is given in 
the Appendix. 



3.2 NGC 2264 G 

The two outflow lobes of NGC 2264 G can be seen just to the 
northeast of S Mon in Fig. [2] The driving source ( [Teixeira et al.| 
|2008[ [Gomez et al.|p^4 . ) has a mass, including the envelope, of 
2-4 Mq. [Teixeira et al. ( |2008| ) report Spitzer imaging of the pro- 
tostellar jet, showing three changes in direction. Fig. [7] shows the 
red and blue contours for three diff'erent velocity ranges through 
the outflow lobes, overlaid on the H2 emission from the protostel- 
lar jet. The CO emission from the molecular outflow follows the 
emission in the jet, including the changes in direction. There is a 
strong correlation between the structure of H2 emission and high 



velocity CO outflow lobes in the NGC 2264 G outflow. The extent 
of the H2 emission matches the size of the CO lobes, and the H2 
intensity peaks are close to the peaks seen in the CO maps. 



The total extent of the red outflow lobe is ~280'', or 1.1 pc 
(assuming 800 pc to NGC 2264), with the very high velocity gas 
(>35 kms~^) extending further than the lower velocity gas. The 
low velocity emission (6 km s~^) traces a cavity, while the high ve- 
locity emission (40 km s~^) traces a collimated flow. The blue lobe 
has a small high velocity extension, and the total extent of the blue 
outflow lobe seen is ~200'', or 0.8 pc. These figures are in agree- 
ment with those derived from CO 2^1 observations ( [Lada & Fich| 
|1996| ). The diff'erence in the highest velocity spatial extent between 
the red and blue lobes we detect in CO 3^2 may be merely the 
amount of material available for entrainment in the two directions. 
The red lobe in Fig. [TJa) suggests a wide opening angle, with the 
H2 emission (MHO 1359) associated with the northern edge of the 
flow. The clumpy structure seen very clearly in our new images of 
MHO 1359 could result from shocks fronts in the boundary layer 
between the flow and the ambient medium. The lack of H2 emis- 
sion along the southern edge of the flow lobe could be simply be 
due to a lower ambient density, suggesting a gradient in the ambient 
density that increases to the north across this region. Alternatively, 
[Teixeira et al.| ( [2008| propose a slowly precessing jet, plus addi- 
tional deflection, to explain the structure. The broad red lobe we 
detect may suggest a density gradient, rather than a precessing jet. 
Table [1] gives the calculated outflow parameters for the NGC 2264 
outflow. The total mass entrained in the outflow is 0.9 M© with a 
momentum of 11 M© kms~^ This is consistent with the mass and 
momentum found for t he lower excitation CO 2^1 transition 
Lada & Fich|l996| ). The jet is bright, and 



(1 M©, 12M©kms- 

has the highest L//2 luminosity for an individual flow in this region. 
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Figure 6. Red and blue shifted CO emission contoured over the H2 continuum-subtracted greyscale image towards the NGC 2264 D cluster. Contours are at 
3.9, 6.5, 9.1, 14.3, 19.5, 27.3, 37.7, 48.1 K km s"^. Red- shifted emission extends from 11.1 to 25.5 kms"\ and the blue-shifted from -18.5 to 1.0 kms"^. 
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Figure 7. Red- and blue-shifted CO emission contours for three different 
velocity ranges towards NGC 2264, overlaid on the H2 greyscale image. 
The red-shifted emission is associated with MHO 1359, the blue-shifted 
emission with MHO 1358. (a) Low-velocity gas at 6.0-8.0 kms~^ (red) 
and 1.0-3.0 kms~^ (blue). CO contours start at 2.5 Kkms~\ with steps of 
2.5 Kkm s"^ (b) Mid- velocity gas at 16.0-26.0 km s"^ (red) and -17.0- -7.0 
(blue), contours as (a), (c) High-velocity gas at 29.0-39.0 km s~^ (red) and 
-30.0- -20.0 km s~^ (blue). CO contours start at 2.0 Kkm s~\ with steps of 
2.0Kkms-^ 



3.3 SMon 

SMon is positioned in front of and travelling towards the 
NGC 2264 cloud I Tauber, Lis & Goldsmith|[T993] ), affecting the 
region with an intense ionizing field. The emission surrounding 



S Mon shows many filaments and arcs, although we do not detect 
any evidence from the presence of extended line wings, or the pres- 
ence of H2 knots, any protostellar outflow activity. The emission 
near SMon is seen across a large velocity range, from 6.7 kms~^ 
to 13.5 kms~\ and is characterized by relatively small linewidths, 
indicating that the complex kinematic structure is not due to proto- 
stellar outflow activity. At the highest red and blue velocities, the 
emission is very clumpy and fragmented. In the ring of emission 
(Fig. [2]), a cluster of SCUBA cores to the north-west is associated 
with CO clumps at blue-shifted velocities, while the SCUBA cores 
to the south-east are associated with CO clumps at red-shifted ve- 
locities. 



3.4 NGC 2264 IRAS-2 

Fig. [8] shows outflow lobes associated with a source first detected 
by Margulis , Lada & Youn g (1989 ). This is a far-IR source with no 
optical counterpart and no previously detected outflow, designated 
as IRAS-2, with a luminosity of 6.7 Lq. |Wolf-Chase, Wlaker &| 
|Lada|(T995] ) detected an extended CS source at this position, but no 
outflow. Along with a weak CO molecular outflow, we also detect 
H2 knots indicating protostellar jet activity, in MHO 3109 (Fig. [8]). 
The outflow candidate NGC 2264 A (Mar gulis, Lada & Snell] 
p^88] ) is -400'' east of IRAS-2. Although a red outflow lobe has 
previously been detected towards NGC 2264 A, we do not clearly 
detect the outflow lobe in CO 3 -^2 emission, or detect any H2 
emission. 



3.5 IRAS 06396+0946 

East of the main emission region, at the position 06^42"^25^5 
+09 43 09 is a cometary- shaped feature, surrounded by several 
small, compact regions in ^^CO, which are blue-shifted to the north- 
west, and red- shifted to the south-east. The peak antenna tempera- 
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Figure 8. Red- and blue-shifted CO emission contours, overlaid on the H2 
greyscale image towards MHO 3109, shown as an ellipse. Red CO contours 
start at 3.0 Kkms~^ with 1.0 Kkms~^ steps, and blue contours start at 2.0 
Kkms~^ with 0.8 Kkms~^ steps. A cross marks the position of the far-IR 
source IRAS2 ( Margulis, Lada & Young|1989) . 
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Figure 9. The isolated cometary feature east of the main NGC 2264 fila- 
ment. Left: ^^CO integrated intensity in greyscale. Right: The velocity field 
shown through the intensity- weighted velocity image. The velocity data has 
been thresholded at the 3cr level for clarity. 



ture is 10.8 K, and the lines are narrow, with FWHM ~1.0 kms~ 
from a Gaussian fit to the data. Fig. [9] shows the integrated in- 
tensity image of this region, in the velocity range 10.1 kms~^ 
to 13.9 kms~^ along with the intensity-weighted velocity image, 
which shows the velocity field. This is an intriguing structure, 
which could be associated with IRAS 06396-^0946, although there 
is as yet no other published data on this region. The velocity field 
suggests an outflow or jet, although the morphology of the ^^CO 
emission is more reminiscent of H2 knots, such as those in the 
blue- shifted jet of NGC 2264 G, rather than the more linear or ellip- 
tical structures more commonly seen in CO. There is no H2 emis- 
sion associated with this object. The spatial alignment and velocity 
structure of this region may just be coincidental, however, with the 
emission arising in areas at diflferent distances which are not re- 
lated. More sensitive observations would be required to determine 
whether these emission regions are connected. 



3.6 The NGC 2264 D filament 

To the east of NGC 2264 D, there are two newly detected MHOs, 
MH01378 and MH01379, shown in the Appendix (Fig.[TT]3). The 
location of the MHOs and the velocity structure of the CO emis- 
sion suggests that these MHOs may be associated with molecu- 
lar outflows which are erupting from the dense gas in the cluster. 



which contains several 5/?/^z^r-identified protostars ( |Sung et al. | 
[2009 1 ). However, the line profiles do not have the line wings usu- 
ally associated with molecular outflows. As the velocity increases 
from 1.0 kms~\ two filaments extend eastwards from the bulk of 
the cluster, eventually overlapping and merging at velocities ~5.0- 
6.5 km s~^ , when the overall length is ~1 T. At higher velocities, the 
filament becomes more pronounced, with additional internal struc- 
ture that is very linear. The brightest regions are those furthest from 
the central cluster. 



4 PRINCIPAL COMPONENT ANALYSIS 

In order to undertake an analysis of the turbulent characteristics 
of NGC 2264, we have carried out a principal component analysis 
(PC A) on the ^^CO data. The method is relatively robust against 
the eflfects of resolution and noise. The low-order components con- 
tain features that contribute most to the variance of the data, while 
higher-order components contain more subtle features within the 
spectral shapes. 

We have implemented the technique described by Heyer &| 
|Schloerb| ( [T997l ); |Brunt & Heyer| i |2002a| l and |Brunt & Heyer| < |2002| | 
who detail the methods in full. This technique transforms the orig- 
inal spectroscopic data cube (RA, Dec and velocity) into a set of 
orthogonal functions which are described by the principal com- 
ponents, /, ordered by decreasing variance of the data projected 
onto each orthogonal vector. The principal components consist of 
an eigenvector, tracing only the velocity structure, and an eigen- 
image, constructed from the projection of the data onto each of 
the eigenvectors. Each eigenimage therefore contains only spatial 
structure, mapping the size scales of diflferences in the line pro- 
files traced by the eigenvectors. The principal components have 
characteristic length scales (L/) determined from the eigenimage, 
and characteristic velocity scales (Svi) determined from the eigen- 
vectors. The eigenvectors and eigenimages of each principal com- 
ponent are coupled, so the physical dynamics can be determined 
from the characteristic scales. The data can be reconstructed from 
a linear combination of only the significant principal components, 
which account for most of the variance, or variability, in the data. 

The technique provides an objective method of extracting the 
most significant components of the data for analysis, the challenge 
then being to interpret the results in a physically meaningful way. 
[Brunt & Heyer^ (2002aJ have carried out PCA on simulated data 
in order to provide empirical relations between the derived PCA 
results and the statistics of the original data, which we utilize in 
this analysis to aid the physical interpretation. 

[Brunt, Heyer & Mac Low (2009) used numerical magneto- 
hydrodynamic models and molecular spectral line observations to 
show that the ratio of the characteristic length scales of the first two 
principal components (L2/L1) is related to the ratio of the turbulent 
driving scale to cloud size, (Ad/Ld) for isotropically forced turbu- 
lence (see also Brunt 2003 ). Therefore, observational measures of 
L2/L1 can be used to estimate Ad/Ld. For ratios > 0.2, there is little 
sensitivity to the actual driving scale, and large-scale driving best 
describes the turbulent driving scale ( [Brunt et al.|2009p . 



4.1 Principal Components 



Following [Brunt[ ( [2003 ), we do not subtract the mean from the data, 
so that /i approximates the integrated intensity and mean line pro- 
file of the ^^CO data. The scale of data values in the eigenvectors 



Wide-field imaging of NGC 2264 9 



I 


Variance 


Cumulative Variance 


Li 


5vi 




% 


% 


pc 


kms~^ 


1 


74.15 


74.15 


7.5 


6.3 


2 


20.86 


95.01 


2.6 


2.5 


3 


3.16 


98.18 


1.3 


1.7 


4 


0.73 


98.90 


0.6 


1.3 


5 


0.51 


99.41 


0.6 


0.8 


6 


0.21 


99.63 






7 


0.08 


99.71 






8 


0.05 


99.76 






9 


0.03 


99.79 






10 


0.02 


99.82 







Table 2. Contribution of the first 10 principal components to the NGC 2264 
CO data. Characteristic length and velocity scales for the 5 significant prin- 
cipal components are listed. 



is not directly related to flux, but to the contribution of the princi- 
pal component to each velocity channel. Large relative positive or 
negative values indicate that a principal component contributes sig- 
nificantly in that velocity channel. The eigenimages can be used as 
a diagnostic of how much spatial structure is associated with each 
eigenvector. 

Table |2]lists the first 10 principal components, along with the 
variance of and the total cumulative variance for each component. 
The first five principal components contribute >99% of the vari- 
ance, with /i accounting for 74% of the variation in the data. The 
largest values of Li corresponds to the largest values of 6vi indicat- 
ing that the largest velocity diflferences are distributed on the largest 
scales. 

Fig. [To| shows the eigenimages and eigenvectors associated 
with the first 5 principal components. The black and white (posi- 
tive and negative) regions in the eigenimages trace velocity fluctu- 
ations whose magnitude is indicated by the width of features in the 
eigenvectors. The first principal component, /i, as noted above, has 
an eigenimage which approximates the ^^CO integrated intensity 
map. The centre velocity of the emission is 7.4 km s~\ and it has a 
large an FWHM linewidth of 6.9 km s~\ from a Gaussian fit to the 
eigenvector, which approximates the mean line profile of the ^^CO 
data. I2 shows a positive-negative dipole pattern that is characteris- 
tic of giant molecular clouds ( [Brunt et al.|2009| ), when the turbu- 
lence is driven on scales comparable to the cloud size. The lower 
velocity contributes mostly to spatial variations in the north, while 
the higher velocity contributes mostly to spatial variations in the 
south. This is consistent with the velocity field we determined for 
NGC 2264, discussed in Sec. |2.3| The positive velocity component 
of I2 is at 5.4 km s~^ , with a linewidth of 4.2 km s~^ . This is spatially 
most closely associated with the emission from the star-forming 
clusters in the cloud, NGC 2264 C and NGC 2264 D. The nega- 
tive velocity component of k is at 10.4 km s~\ and has a narrower 
linewidth of 2.6 kms~^ This component is spatially most associ- 
ated with emission from the region near S Mon, and the tip of the 
Cone, /s, with peaks at 4.0, 8.2 and 10.9 kms~\ and linewidths of 
2.9, 1.7 and 1.6 kms~\ describes velocity and linewidth variations 
within the velocity components identified in /i and I2. All of the 
significant principal components also show an isolated, compact 
region of emission to the east of the main emission region, at the 
position associated with IRAS 06396+0946 (discussed in Sec. [3]). 

Further principal components describe successively smaller 
spatial variations in velocity and linewidths. By /20, there is negli- 
gible emission in the eigenimages, and from these high-order prin- 



cipal components, noise estimates can be made ( [Brunt & Heyerj 
|2502t . 



4.2 Characteristic Scales 

Using results from the PCA decomposition, we can investigate the 
characteristic velocity variations and the spatial scales over which 
they occur. The characteristic scales are calculated using the nor- 
malized autocorrelation functions (ACFs) of each principal com- 
ponent, /, following [Brunt & Heyer[ ( [2002a[ ). The spatial scales 
are calculated using noise- subtracted, resolution-corrected ACFs. 
Given the limited number of spectral pixels, the characteristic ve- 
locity scale measurements are largely unaflfected by instrumental 
noise. Instrumental noise and finite resolution eflfects are explained 
in detail by Brunt & Heyer[ ( [2002a[ ), whose procedures we have im- 
plemented. The characteristic velocity and length scales (6vi and 
Li) are determined from the lag at which the normalized ACFs of 
the eigenvectors and eigenimages fall to ^~\ and are calculated fol- 
lowing [Heyer^Schloe^ ([199^ 

The characteristic length scales L/, and the characteristic ve- 
locity scales Svi for the 5 significant principal components are listed 
in Table [2] Both the length and velocity scales decrease with in- 
creasing /. These scales can be visualized in the eigenimages and 
eigenvectors, where we can see, respectively, the spatial scale of 
variation, and the width of spectral features, decreasing as / in- 
creases. The assumption of a PCA analysis of molecular cloud 
spectral line data is that the decreasing size scales are reflecting a 
fundamental property of the kinematics of the gas within the cloud 
( [Heyer & Schloerb[1997| l, from which turbulent driving scales (Ad) 
can be determined. The longest flow we detect towards NGC 2264 
is that of NGC 2264 G, at 280'', or 1.1 pc. Comparing this value 
with the L/ size scales suggests that the outflows we detect are con- 
tributing to the variance in the data at characteristic scales < L3, and 
so are contributing <3% to the variance in the data. They are not, 
therefore, contributing at the size scales necessary to be the main 
component driving turbulence. This is supported by our analysis in 
Sec. [3] where we found that <0.5% of the cluster masses have been 
entrained within molecular outflows. 

On the largest scales, the velocity variations are due to dif- 
ferences between components that make up the entire molecu- 
lar cloud. The characteristic spatial scale of the first eigenimage, 
Li-1 .5 pc, therefore gives an estimate of the overall cloud size, Ld. 
The characteristic scale of the second eigenimage, L2-2.6 pc, mea- 
sures the size- scale over which the largest velocity variations occur 
in the data. L2-2.6 pc can therefore be associated with the turbulent 
driving scale. 

[Carroll, Frank & Blackman] ( [2010[ ) have examined the use 
of PCA in simulations of outflow driven turbulence, isotropically 
forced turbulence, and a combination of both. They find that the 
ratio of characteristic length scales /2//1 correctly identifies the 
largest driving scale in all three types of simulations, but suggest 
the scales derived from PCA are measuring the largest scale of co- 
herent motion, which is not necessarily the turbulent driving scale. 
As noted by [Carroll et al.[ ( [2010[ ), this distinction is important in 
clouds where outflows are present within an external cascade, or 
where there are large-scale coherent motions due to angular mo- 
mentum conserving collapse. In our PCA analysis, the principal 
components describe a large-scale velocity field across the cloud, 
with contributions from high velocity outflow emission. As shown 
in Fig. [4] the velocity field of NGC 2264 shows little evidence for 
large-scale coherent motions. There is no indication of a rotation 
axis for the region as a whole. 
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We follow [Brunt et al.| { |2009| ) and use the ratio of the charac- 
teristics length scales, /2//1, to estimate the turbulent driving scale. 
From the ratio of the spatial scales of the first two eigenimages, we 
derive a fractional driving scale of 0.35, indicating large-scale driv- 
ing. Combined with the characteristic scale L2-2.6 pc associated 
with the turbulent driving scale, and our analysis in Sec. |3] these 
data indicate that protostellar outflow activity in NGC 2264 is not 
the dominant component of turbulence, even though it is a young 
and energetic star forming region. 



4.3 Size-linewidth correlation 



|Larson| ( |1981| ) developed an empirical relation between the 
linewidth AV and cloud size R, with AV oc R", and a-0.5. This 
power-law relation is similar to that expected from studies of tur- 
bulence, where values of a can range from of = 1/3 for the in- 
compressible energy cascade io a = 1 12 for a shock-dominated 
turbulent fluid ( .McKee & Ostriker|2007| and references therein). 
The equivalent PCA velocity statistic a is obtained from the set of 
significant 6vi - Li pairs which are larger than the resolution limits. 
These form the power-law relationship dv oc L"^. [Brunt] (2003) com- 
pared observations to simulations, and obtained values of a from 
0.5 t o 0.8. 

[Brunt et aT] ( [2003| ) relate the exponent of the characteristic 
velocity and length scale relation, a, to the scaling properties of 
the intrinsic velocity field. Roman-Duval et al.|([2011) derive a 
PCA calibration for the exponent of the turbulent velocity spectrum 
E{k) oc r^^ finding /5v = 0.2 + 0.05 + (2.99 + 0.09)Qf. Values of 
)S = 5/3 are expected for the incompressible energy cascade model 
([Elmegreen & Scale|2007| and references therein). [Roman-Duval et[ 
|al.[(|201 1|) applied their results to a sample of 367 molecular clouds, 
finding {a) = 0.61 + 0.2, and (/?> = 2.06 + 0.6. 

We calculate a = 0.74 + 0.08 and /5 = 2.4 from our 
PCA analysis. The value of a is larger than that expected from a 
classic linewidth- size relationship (Larson 1981) in giant molec- 
ular clouds, but is consistent with measurements towards the 
NGC 2264 C region ( [Maury et al.|[2009] ) using observations of 
^^CO. The size-linewidth correlation is thought to depend upon the 
mean surface density of the cloud ( jHeyer et al.[|2009| ), and so a 
values will be higher in regions of high surface density. NGC 2264 
has a mean surface density >10 times that assumed for the classic 
linewidth- size relation ( [Maury et al.|2009| ). a - 0.74 is consistent 
with values found in high mass star forming regions such as the 
Rosette molecular cloud ( jHeyer, Williams & Brunt 2006). 

The p value of 2.4, giving E{k) oc k~^-^, is larger than values 
expected for Kolmogorov-type turb ulence, but is sim ilar to values 
found using ^^CO towards Perseus ( Sun et al.|2006 l, where values 
were higher in active star forming regions than in dark clouds. A 
steep energy spectrum is consistent with the energetic star forma- 
tion activity of NGC 2264. 



5 SUMMARY 

We have presented wide-field spectral imaging observations in 
^^CO 3^2 and wide-field high resolution imaging observations 
in H2 1-0 S(l) towards NGC 2264. These observations, covering 
nearly 1 square degree, off'er a detailed view of the star formation 
activity taking place in this clustered environment. We find that 

(i) Protostellar outflow activity in NGC 2264 does not occur on 



large enough scales to drive the turbulence. The largest flows ex- 
tend to ~1 pc, while a PCA analysis suggests that turbulence is 
driven on scales larger than 2.6 pc. 

(ii) Only a small fraction, <0.5%, of the cluster gas mass in 
NGC 2264 C and NGC 2264 D, has been swept up to high veloci- 
ties through protostellar outflow activity. The outflow activity is not 
having a sufficient impact on the cloud to be the dominant source 
of turbulence. 

(iii) We detect 46 molecular jets and knots in H2, 35 of which 
are new detections. Based on the new H2 data, NGC 2264 C and 
NGC 2264 D extend to a larger spatial extent than previously deter- 
mined. 

(iv) NGC 2264 D contains a larger mass of gas at higher veloci- 
ties, and also more energy and momentum than NGC 2264 C. 

(v) Due to the presence of spatially distinct red- and blue-shifted 
CO emission, and H2 emission in the form of jets or bow shocks, 
we characterize 4 cores in NGC 2264 C as protostellar that were 
previously identified as prestellar. 

(vi) Of three massive YSO candidates in the NGC 2264 D clus- 
ter, we detect a bipolar molecular outflow and H2 jet towards 
SSB 11829. 

The large number of protostars, jets and outflows in the 
NGC 2264 region, and their small separations, require very high 
resolution observations in order to isolate the flows and identify 
driving sources. This will be possible in future observational stud- 
ies with the resolution and sensitivity off'ered by ALMA. 
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7 APPENDIX A 

7.1 Molecular Hydrogen emission-line Objects 

Table [3] lists the 46 Molecular Hydrogen emission-line Objects 
(MHOs) identified in this work. Of these, MHO 1375-1399 and 
MHO 3100-3109 are new discoveries. Fig. [TT|fT2l show the H2 
emission maps of the MHOs listed in Tablel3] H2 emission maps of 
NGC 2264 C and NGC 2264 D are shown in the next section to aid 
visualization when discussing the clusters. 

7.2 The protostellar content of NGC 2264 C 

NGC 2264 C is associated with IRSl, or Allen's source, an early B 
star discovered by |Allen| ( |1972| ), who estimated extinction towards 
the source as A^ ~26. The group of prestellar cores and protostars 
to the south are known as the cluster NGC 2264 C ( Margulis, Lada| 
|& Snellj 1988) [Peretto et al.|2006| ), where star formation has prob- 
ably been triggered by activity from IRS 1 ( [Thompson et al.|1998[ 
[Williams & Garland||2002|l. The cores are forming several inter- 
mediate/high mass stars (Peret to et al.] [2006, Ward-Tho mpson et[ 
[^..200 0), with core masses of 1-50 M© measured from submil- 
limetre observations. Previous studies with a variety of molecular 
lines have pointed to the complex velocity field in this cluster, with 
IRSl associated with velocities between 8-10 km s~ \ surrounded 
by several sub-clouds with velocities of 5.5-7 kms~^ |Schreyer et 



|ai[1997j . The emission also extends to the west, with no obvious 
driving source, or any H2 counterparts. IRSl itself drives a tightly 
peaked, overlapping molecular outflow. There are several submm 
and Class I cores clustered in the region of the IRSl source that 
could be driving these energetic flows. 

[Peretto et al.[ p006| l identifies 4 prestellar cores, and 8 pro- 
tostellar cores towards the clump near NGC 2264 C, from ob- 
servations of the 1 .2 mm dust continuum, and the optically thin 
tracer N2H+. Fig. [13] shows the CO molecular outflows and H2 
jets, marked with the positions of the pre- and protostellar cores, 
adopting the labelling scheme of [Peretto et aL[ ( [2006[[2007| ), and 
labelling the cores CMM1-CMM13. Of these, we detect possible 
outflows towards all of the prestellar cores, based on the presence of 
spatially distinct red- and blue-shifted emission, and H2 emission 
in the form of jets or bow-shocks. The presence of H2 jets indi- 
cates that these cores are young Class sources. The cores CMM6, 
CMM7 and CMM8 are located to the south of NGC 2264 C, and 
[Peretto et al.[ ( |2006| ) indicates that they are not associated with 
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Table 3. Molecular Hydrogen emission line Objects (MHOs) in NGC 2264*. 



Object'^ 



RA 

(2000.0) 



Dec 

(2000.0) 



Assoc.^ 
HH Object 



Note 



MHO 1349 


06 41 10.5 


+09 29 47 


MHO 1350 


06 41 09.9 


+09 28 14 


MHO 1351 


06 41 12.2 


+09 29 50 


MHO 1352 


06 41 15.3 


+09 29 44 


MHO 1353 


06 41 17.9 


+09 30 08 


MHO 1354 


06 41 15.3 


+09 28 34 


MHO 1355 


06 41 17.7 


+09 27 59 


MHO 1356 


06 41 12.7 


+09 28 03 


MHO 1357 


06 41 16.8 


+09 27 10 


MHO 1358 


06 4105.0 


+09 55 50 


MHO 1359 


06 4123.0 


+09 55 45 


MHO 1375 


06 4104.2 


+09 46 35 HH 125 


MHO 1376 


06 4103.7 


+09 44 32 HH 225 


MHO 1377 


06 4102.2 


+09 39 51 HH226 


MHO 1378 


06 4136.5 


+09 36 31 


MHO 1379 


06 4129.3 


+09 36 08 


MHO 1380 


06 41 14.3 


+09 37 50 


MHO 1381 


06 41 10.6 


+09 38 10 


MHO 1382 


06 41 11.3 


+09 36 42 


MHO 1383 


06 4104.8 


+09 36 07 


MHO 1384 


06 4102.7 


+09 35 55 


MHO 1385 


06 4100.0 


+09 36 19 


MHO 1386 


06 41 14.3 


+09 35 35 


MHO 1387 


06 41 14.2 


+09 34 10 


MHO 1388 


06 4105.6 


+09 34 02 


MHO 1389 


06 4104.8 


+09 34 44 


MHO 1390 


06 4106.2 


+09 33 44 


MHO 1391 


06 4104.2 


+09 33 32 


MHO 1392 


06 40 59.4 


+09 33 50 


MHO 1393 


06 40 54.9 


+09 33 55 


MHO 1394 


06 40 51.9 


+09 34 18 


MHO 1395 


06 40 52.9 


+09 33 19 


MHO 1396 


06 4101.5 


+09 3157 HH581 


MHO 1397 


06 40 48.2 


+09 32 58 


MHO 1398 


06 40 47.2 


+09 32 13 


MHO 1399 


06 40 39.8 


+09 3141 


MHO 3100 


06 40 33.6 


+09 32 00 


MHO 3101 


06 40 52.3 


+09 29 43 


MHO 3102 


06 4101.4 


+09 30 56 


MHO 3103 


06 41 10.2 


+09 30 27 


MHO 3104 


06 41 19.9 


+09 29 41 


MHO 3105 


06 4128.8 


+09 28 24 


MHO 3106 


06 4107.4 


+09 30 04 


MHO 3107 


06 4106.8 


+09 26 54 


MHO 3108 


06 4109.6 


+09 25 45 


MHO 3109 


06 42 09.5 


+09 13 09 



Emission associated with NGC 2264 C 

Knotty jet extending to south of NGC 2264 C/eclipsing TTS KH 15D 

Arc in NGC 2264 C region; possibly associated with MHO 1350 

CoUimated jet 1-2 arcmin ESE of NGC 2264 C 

Chain of knots in the NGC 2264 C region 

Compact knot in the NGC 2264 C region 

Small knot in the NGC 2264 C region 

Compact feature in the NGC 2264 C region 

Small, faint knot in the NGC 2264 C region 

Complex group of knots associated with bipolar outflow NGC 2264 G 

Faint arcs and filaments associated with bipolar outflow NGC 2264 G 

Knots and extended bows ~7' south of S Mon in NGC 2264 

Knots and faint emission ~7' south of MHO 1375 (unrelated) 

Two compact knots to the north of a cluster of outflows near NGC 2264 D 

Compact knot to the west of the NGC 2264 D cluster 

Faint feature to the west of the NGC 2264 D cluster 

Faint knots and emission in NE corner of the NGC 2264 D cluster 

Compact knot plus difl'use emission in NE corner of the NGC 2264 D cluster 

Bright, knotty arc plus difl'use emission in the NGC 2264 D cluster 

Faint emission in the NGC 2264 D cluster 

Faint emission in the NGC 2264 D cluster 

Bright knot and extended west-facing bow shock in 

Bright, knot bows moving north-eastward from the NGC 2264 D region 

Arcs and filaments in flow moving eastward from the NGC 2264 D region 

Compact knots and fingers off' emission in the NGC 2264 D cluster 

Knots and emission in the NGC 2264 D region 

Extended, collimated outffow in the NGC 2264 D cluster 

Bright, extended outffow that crosses MHO 1390 in the NGC 2264 D region 

Faint ffngers of emission extending to the south-west in NGC 2264 D 

Group of knots and arcs in the NGC 2264 D region 

Smaff group of knots; possibly associated with the MHO 1395 jet 

Bright bow shock and knotty jet to the west of the NGC 2264 D cluster 

Very faint outffow in the NGC 2264 D region 

Knot and faint emission a few arcmins to west of the NGC 2264 D cluster 

Knot and faint emission a few arcmins to west of the NGC 2264 D region 

Arcs plus diff'use emission in feature to west of the NGC 2264 D region 

Small group of knots plus faint emission -1.5' west of MHO 1399 

Compact feature to south of the NGC 2264 D cluster 

Compact feature to south-west of the NGC 2264 D cluster 

Curving chain of knots/outff ow ~1' north of NGC 2264 C 

Knot plus faint ffngers of emission ~2.5' east of NGC 2264 C 

Diff'use arcs of emission -4.5' ESE of NGC 2264 C 

Two compact knots ~V north-east of NGC 2264 C 

Knot plus faint emission ~3' south of NGC 2264 C 

Bright knot and arc of emission ~4' south of NGC 2264 C 

Chain of half-a-dozen knots in a ffow in the south of NGC 2264 



^ For discovery images of MHO 1349-1357 see Wang et al. (2002); MHO 1358/1359 was ffrst imaged in H2 by Davis & Eisloff'el (1995). MHO 1375-1399 

and MHO 3100-3109 are new discoveries. 

^ Associated HH object, if any (see Reipurth et al. 2004, for details). 

* For a complete list and individual images of MHOs, see: ht tp://www.jach.hawaii.edu/UKIRT/MHCat/| 



any MSX or 2MASS sources, or any H2 jets. In Fig. [13] CMM7 
is associated with a relatively bright Spitzer 24 \im source, while 
CMM6 is spatially coincident with a faint 24 \xm source, and all 
three are associated with red- and blue- shifted elongated emission 
which could be arising in protostellar outflows. In addition, CMM9, 
the northern-most source to the east of NGC 2264 C also appears 



to have a molecular outflow, and be associated with a weak 24 \im 
source. We propose that these are all, in fact, protostellar sources, 
induced by the activity of NGC 2264 C. Table. |4] summarizes the 
CO outflows and MHOs we have been able to associate with a driv- 
ing source, and the lengths of the flows, measured from either the 
H2 emission, or from CO red- or blue-shifted emission. 
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Figure 11. H2 colour composite images (as Fig. [3) of a) NGC 2264 G. b) The NGC 2264 D filament, c) NGC 2264 IRAS2 
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Figure 12. H2 colour composite image of S Mon, as Fig. [3] 



CMMl drives a double- sided H2 jet, MH01352, and has a 
clearly defined red outflow lobe. The blue outflow lobe is driving 
towards the centre of the cluster, directly towards CMM9, and so 
can only be seen as a spatial extension in the blue shifted emission, 
aligned along the H2 jet direction. 

CMM2 is at the base of a chain of H2 knots and fingers. 



MHO 1304, that is also associated with weak blue-shifted emis- 
sion. We are not able to associate any red-shifted emission from 
this source due to the strong outflow lobes at the centre of the clus- 
ter. A Spitzer 24 |^m identified protostar, 14214 ( |Sung et al. |2009| l 
is8''northofCMM2. 

MHO 1356 is associated with a large red- shifted outflow lobe. 
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Figure 13. Top: H2 colour composite image of NGC 2264 C, as Fig. [3] Bottom: Red and blue shifted CO emission contoured over the H2 continuum- subtracted 
greyscale image. Contours are at 3.9, 6.5, 9.1, 14.3, 19.5, 27.3, 37.7, 48.1 Kkms~^ Red- shifte d emission e xtends from 12.2 to 30.9 kms~\ and the blue- 
shifted from -17.8 to 2.1 kms~^. The saturated source seen in the H2 emission is Allen's source JAllen 19721. The sources identified by 1 
are marked and labeled with the CMM number. 



Peretto et al. 



2006t 



The jet, shows a curved chain of features, extends southwards from 
CMM3, and so we tentatively assign this as the driving source. The 
lobes can't be traced directly back to the source due to the overlap- 
ping flows in the central regions. 

CMM7, previously identified as prestellar, is at the base of a 
bipolar molecular outflow and associated bow shocks, MHOS 108 
and MHO 1354. CMM7 is also associated with a Spitzer 24 |j,m 
identified protostar, 13822 ( |Sung et al."||2009l ). CMM6 is co- 



incident with a chain of H2 knots ending in a bow shock in 
MHO 1355, associated with red-shifted emission. There is also a 
small extension of H2 knots seen associated with blue-shifted emis- 
sion on the other side of the source. If this jet extends any fur- 
ther, the emission is not distinguishable from the emission extend- 
ing from the Cone Nebula. CMM8 is located between two small 
red- and blue-shifted emission regions which it may be driving, al- 
though there is also an H2 emission knot perpendicular to the line 
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Source 




H2 


CO lobes 


length {") 




Red 


Blue Red 


Blue 


Red 


Blue 


CMMl 


MHO 1352 


Y 


Y 


45 


58 


CMM2 




MHO1304 ? 


Y 




92 


CMM3 


MH01356 




Y 


? 


75 




CMM6 


MH01355 




Y 


Y 


67 


14 


CMM7 


MHO3108 


MHO 1354 Y 


Y 


90 


100 


CMM8 


9 




? Y 


Y 


31 


20 


CMM9 


MH01353 




Y 


? 


77 




CMM10/13288^MHO3106 




Y 


? 






CMM13 


MH01351 


MHO1350 Y 


Y 


126 


124 


13350^ 


MHO 1303 


Y 


N 


41 


46 


14314714415^ 


MH01357 


7 


7 






^ Identification from 


Sung 


etal. 


2009 











Table 4. Driving source identifications for NGC 2264 C flows. 



of the outflow lobes, at a similar distance to the source. There is 
a Spitzer-i&Qn\ifiQ& protostar at this location, 14314. MHO 1357 is 
spatially co-incident with the two Class I Spitzer identified sources, 
14314 and 14415 jSung et al. |[2009] ) separated by just lr^ The 
sources and outflows in this region are closely clustered, and driv- 
ing jets and outflows in dififerent directions. This small cluster of 
cores (CMM6, CMM7 and CMM8) is surrounded by compact H2 
emission, and it is difficult to definitively assign all of the emission 
to one source or another. The amount of H2 emission and red- and 
blue-shifted CO emission does suggest that these cores are proto- 
stellar, rather than prestellar as previously catalogued. 

CMM9 has been identified as prestellar by ( [Peretto et al.| 
|2006| l, but we see evidence of a chain of H2 knots, MHO 1353, 
and weak red-shifted emission associated with this source, perpen- 
dicular to the blue- shifted outflow lobe from CMMl. Any blue- 
shifted emission that may be associated with this source is entan- 
gled with the blue outflow lobe from CMMl. There is a Spitzer 
24 i^m identified protostar, 14214 < |Sung et al."||2009l ) IV from 
CMM9, which is spatially coincident with the terminal bow shock 
from the blueOshifted lobe of CMMl . 

CMMIO is close to a Spitzer 24 |j,m Class I source (13288 
[Sung et al. ||2009| ), very close to IRSl, which is saturating the 
H2 image. The Spitzer peak is spatially most closely aligned with 
the red outflow lobe, and is close, although not coincident with 
CMMIO. There is also two compact H2 knots, MHO1306, which 
align with this source and red-shifted emission. Therefore, we ten- 
tatively identify CMMIO as the driving source of this flow. Any 
blue-shifted emission is co-incident with the energetic flows at the 
centre of the cluster, and not easily distinguishable. 

CMMl 3 has a spectacular bipolar jet and molecular out- 
flow, and has a particularly clear jet and multiple bow- shocks, 
MHO3107, associated with the blue-shifted lobe. The red-shifted 
lobe is coincident with MH01351. 

The remaining CMM sources are all tightly clustered near 
Allen's Source, which drives an energetic bipolar flow, and is satu- 
rated in the H2 emission, making it difficult to establish individual 
flows. 



7.3 The protostellar content of NGC 2264 D 

NGC 2264 D is less luminous than NGC 2264 C, forming several 
intermediate/high mass stars ( [Peretto et al.|2006l ) with core masses 
of 1.9-17.3 Mq, estimated from submillimetre observations. The 
submillimetre sources are less clustered than in NGC 2264 C, but 



high resolution observations ( [Sung et al. |2009[[Teixeira et al.|2006| ) 
show many of the submillimetre sources are co-incident with mul- 
tiple tightly clustered protostars. |Forbrich et al.| ( |20l6 ^) detects PAH 
emission surrounding 3 embedded YSOs in NGC 2264 D, making 
them candidate massive YSOs. We investigate the protostellar sta- 
tus of these objects using our CO and H2 data. All of the candidates, 
SSB 1070, SSB 11829 and SSB 12820 are associated with compact 
K-band continuum sources. SSB 11829 is co-incident with a ^^CO 
bipolar molecular outflow, and MHO 1385, a bright knot and ex- 
tended west-facing bow shock, indicating the source is protostellar. 
A second, SSB 10710, has a very faint H2 emission knot to the west, 
although this is not a sufficiently significant detection to be identi- 
fied as an MHO, and we do not detect any clear outflow signatures 
in ^^CO that we can identify with this source. 

Eleven prestellar and four protostellar cores have been identi- 
fied by [Peretto et al.| ( [2006 '), and we have adopted their labelling 
scheme of, labelling the cores DMM1-DMM15. Fig. [14] shows 
the red- and blue shifted contours overlaid on the H2 emission, 
as Fig. [6] with the DMM sources also labelled. Six of the prestel- 
lar cores are co-incident with Spitzer 24 \im identified protostars 
( [Sung et al. ||2009| ), and two of these show clear outflow struc- 
ture in CO emission, DMM2 and DMM13. Both of these are part 
of a ridge containing a previously identified protostellar source, 
DMM7. The final source in this ridge, DMM 15 is also associ- 
ated with a 24 \im source, but only has weak evidence of outflow 
activity. The NGC 2264 D cluster is more tightly clustered than 
NGC 2264 C, and many of the DMM sources are co-located with 
multiple Spitzer identified protostars, making it diflficult to isolate 
individual flows. Where this has been possible, we list in Table. |5] 
the CO outflows and MHOs we have been able to associate with 
a driving source, and the lengths of the flows, measured from ei- 
ther the H2 emission, or from CO red- or blue-shifted emission. We 
describe the flows in more detail below. 

The area surrounding MHO 1385 contains two prestellar cores, 
DMM4 and DMM6, and 4 5/7/^z^r-identified protostars ( Sung et al. | 
|2009| ), one of which is co-incident with an H2 knot. MHO 13 85 is 
associated with mainly blue-shifted emission, with a more com- 
pact, slightly overlapping red- shifted CO flow to the west. Al- 
though there are multiple sources in this region, the bipolar flow is 
not obviously associated with any of the nearby sources or H2 emis- 
sion. Midway between DMM4 and DMM6 lies a candidate massive 
YSO, SSB 11829 ( [Smith et al.|20T0t , and it is this source that is 
spatially aligned with the molecular outflow and MHO 1385. 50'' 
to the south-east of MH01385 are two more H2 objects, MH01383 
and MH01384, two prestellar cores, DMM2 and DMM15, and 5 
Spitze r-ideniified protostars, all within a radius of 30''. There is 
red- and blue-shifted emission associated with this region, but it 
is not possible to disentangle individual flows at this resolution. 
DMM7 is located at the base of MHO 1386, which is clearly as- 
sociated with blue-shifted emission, although we can't isolate any 
red- shifted counterpart to this flow. 

MHO 13 80 is coincident with a compact peak of red-shifted 
CO emission. MHO 1381 is associated with weak blue-shifted 
emission, and could be a counterpart. There is no known source as- 
sociated with this flow, nor with MHO 13 82, which ends in a bright, 
knotty arc of emission. The shape of the arc suggests it is driving 
towards the cluster, not away from it. MHO 13 82 is associated with 
an extended red- shifted flow, with no detectable driving source. 
MHO 1391 contains multiple arcs and knots of H2 emission, each 
of the H2 clusters associated with a peak of red-shifted CO emis- 
sion - four in total. At the southern end of MHO 1391, there is a 
compact blue-shifted emission region that overlaps the penultimate 



Wide-field imaging of NGC 2264 17 



40'00" 



^VIHtslBSl' 



o 
o 

o 

c 
o 

c 
"u 
Q 



38 '00" 



36 '00" 



34'00" 



""°^¥ (%H0i3a2 



MHQA385 



31386.' ^ *MhI0138B - , - MHOlSg/ 



^MlKli387 :-^ 



JiT^ 



_«U .■ 



MH0139^ 



+9^32^00" 



/ . :.^" - ' . ' MHoi3g9 > 



^pHH03igO 



205 



lOs 41m005 50s 

Right ascension (J2000) 



40s 6h40m30s 



o 
o 
o 
rvi 

c 
o 

"■4-1 

_c 
T3 

QJ 

Q 



40'00" 



38'00" 



36'00" 



34'00" 



+ 9°32'00"- 




20s 



10s 41m00s 50s 

Right ascension (J2000) 



40s 



6h40m30s 



Figure 14. Top: H2 colour composite of NGC 2264 D, as Fig. |3] Bottom: Red and blue shifted CO emission contoured over the H2 continuum- subtracted 
greyscale image. Contours are at 3.9, 6.5, 9.1, 14.3, 19.5, 27.3, 37 .7, 48. 1 Kkm s~^ Red-shifted emission extends from 1 1 . 1 to 25.5 km s~\ and the blue-shifted 
from -18.5 to 1 .0 km s~^ . The sources identified by Peretto et al.|J2006| are marked and labeled with the DMM number. 



red-shifted region. However, there is no source that is obviously 
driving these flows. There are no detected sources at the southern 
end of MHO 1391, although there is a prestellar core, DMM8 to 
the north, and seven Spitzer-idQntifiQd protostars to the north-east, 
near DMMl. DMMl, located within 8'' of 5 5/?/^z^r-identified pro- 
tostars, is also coincident with MH01388. 

MHO 1392 is part of the dominant region of red- shifted emis- 



sion that extends across ~8' to the south-west of NGC 2264 D. It 
is crossed by a double-peaked region of blue-shifted emission, al- 
though there is no nearby source detected that could be driving this 
flow. MH01393, MH01395 and MHO 1394 are also coincident 
with this long ridge of red-shifted emission, although the different 
alignments of the knots, jets and bow shocks associated with these 
MHOs indicate that they may by arising in smaller flows that we are 
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Source 


H2 


CO lobes 


leng 


th n 


Red 


Blue 


Red 


Blue 


Red 


Blue 


DMM7 


MH01386 


? 


Y 




45 


MHO1380 


MH01381 


Y 


Y 


42 


17 


MH01382 




Y 




66 




MHO3100 




Y 


Y 


35 


52 


12253^ 


MHO 1377 


M 


Y 




19 


^ Identification from 


Sung et al. 


2009 









Table 5. Driving source identifications for NGC 2264 D flows. 



not able to distinguish from the larger scale flow. Further along the 
large red- shifted flow, are several more MHOs, MHO 1397, MHO 
1398 and MHO 1399, each of which are coincident with peaks in 
the red- shifted flow. At the south-west tail of the large red flow 
is a compact bipolar flow, that has MHO3100 associated with the 
red-shifted lobe. There is no detected driving source for this flow, 
but it does lie outside the SCUBA and Spitzer observed areas. To 
the north of NGC 2264 D, MH01377 is coincident with a weak 
compact clump of blue-shifted emission, and a Spitzer-i&Qn\ifiQ& 
protostar, 12253 ( |Sung et al. ,2009 ) . 



